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ABSTRACT 



The Jovian dust streams are high-speed bursts of 
submicron-sized particles traveling in the same direction 
from a source in the Jovian system. Since their discov- 
ery in 1992, they have been observed by three spacecraft: 
Ulysses, Galileo and Cassini. The source of the Jovian 
dust streams is dust from Io's volcanoes. The charged and 
traveling dust stream particles have particular signatures 
in frequency space and in real space. The frequency- 
transformed Galileo dust stream measurements show dif- 
ferent signatures, varying orbit-to-orbit during Galileo's 
first 29 orbits around Jupiter. Time-frequency analysis 
demonstrates that Io is a localized source of charged dust 
particles. Aspects of the particles' dynamics can be seen 
in the December 2000 joint Galileo-Cassini dust stream 
measurements. To match the travel times, the smallest 
dust particles could have the following range of param- 
eters: radius: 6 nm, density: 1.35-1.75 g/cm 3 , sulfur 
charging conditions, which produce dust stream speeds: 
220\450kms~ 1 (Galileo\Cassini) and charge potentials: 
5.5\6.3 V (Galileo\Cassini). 

Key words: Jovian dust streams, frequency analysis, 
Cassini-Galileo dust measurements, Io's volcanoes. 



on its way to Saturn, recorded more high-speed streams 
of submicron-sized particles from the Jovian system. The 
many years-long successful Jovian dust streams observa- 
tions reached a pinnacle on December 30, 2000, when 
both the Cassini and Galileo dust detectors accomplished 
a coordinated set of measurements of the Jovian dust 
streams inside and outside of Jupiter's magnetosphere. 

Indirect methods applied by previous researchers have 
pointed to Io being the simplest explanation for the ques- 
tion of the origin of the Jovian dust streams. We first 
show by direct methods that Io is the source of the Jo- 
vian dust streams. To address the issue of identifying 
Io directly in the Galileo dust detector data, we apply 
time-frequency analysis, in particular, Fourier methods, 
to the Galileo dust data. Additional frequency signa- 
tures, such as amplitude modulation, also emerge from 
the time-frequency analysis. 

The second part of this paper focuses on the dust streams 
dynamics. Here, we apply a detailed Jovian particles and 
fields model to simulate a dust stream particle's trajectory 
as the particle moves from Io's orbit through Jupiter's 
magnetosphere and beyond. Through the model, we 
show one possible set of parameters that match the travel 
times seen in the December 30, 2000 Galileo-Cassini 
joint dust stream measurements. 



1. OVERVIEW 



2. IO'S FREQUENCY FINGERPRINT 



The Jovian dust streams are high-speed collimated 
streams of submicron-sized particles traveling in the 
same direction from a source in the Jovian system. They 
were discovered in March 1992 by the cosmic dust detec- 
tor instrument onboard the Ulysses spacecraft, when the 
spacecraft was just past its closest approach to Jupiter. 
Observations of the Jovian dust stream phenomena con- 
tinued in the next nine years. A second spacecraft, 
Galileo, now in orbit around Jupiter, is equipped with 
an identical dust detector instrument to Ulysses' dust in- 
strument. Before and since the Galileo spacecraft's ar- 
rival in the Jupiter system in December 1995, investi- 
gators recorded more dust stream observations. In July 
and August 2000, a third spacecraft with a dust detector 
(combined with a chemical analyzer), Cassini, traveling 



In order to find Io's frequency 'fingerprint' in the Galileo 
dust detector data, we followed the following steps: 1 ) we 
transformed the Galileo dust detector data into frequency 
space via periodograms, 2) we noted frequency patterns 
such as amplitude modulations, 3) we compared the 
frequency-transformed Galileo data with synthetic data, 
and, 4) we noted spacecraft effects such as Doppler shifts. 



2.1. Time-Frequency Analysis via Periodograms 



The "classic" or Schuster periodogram ( Bretthorst 1988 ) 
is conventionally defined as the modulus-squared of the 
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discrete Fourier transform. If the input time series con- 
tains a periodic feature, then the periodogram can be cal- 
culated for any frequency and it displays the presence of a 
sinusoid near one frequency value as a distinct peak in the 
spectrum. The Lomb-Scargle periodogram applied here 
is a slightly modified version of the classic period ogram 
giving a simpler statistical behavior (Scargle 1982). 



Our best Galileo dust dataset for detecting Io's frequency 
fingerprint emerged from the earlier Galileo orbits around 
Jupiter, because the spacecraft orbital geometry of the 
first years of the Galileo mission favored higher fluxes 
of dust stream particles. Therefore, by combining two 
years of the early data, we gained a higher signal-to-noise 
dataset. In Fig. y], we show a Lomb-Scargle periodogram 
for the first two years, 1996-1997, of Galileo dust impact 
rate data. This particular periodogram is the best example 
from the Galileo dust detector data showing, with high 
confidence: Io's frequency signature, Jupiter's frequency 
signatures, and amplitude modulation effects. The peri- 
odogram shows the following frequency signatures. 



Frequency Summary of Fig. [j] 



then the spectrum in frequency space would appear like 
the spectrum in Fig. [j], where the modulation products 
(sidelobes) at Jupiter's frequency at full and half -rotations 
are due to Jupiter's frequency modulating Io's frequency 
of orbital rotation. The frequency difference between 
Jupiter's rotational frequency and each of the sidelobes 
is the same frequency as Io's frequency of orbital rota- 
tion. In addition, if one or both of the original signals 
have broad spectra, then the spectrum of the modulation 
products will be broadened to the same extent (since the 
modulation process is linear). Therefore, the spread of 
Io's peak is repeated in the same way for the sidelobes of 
Jupiter's frequency peaks (full and half-rotation). 

Modulation products appeared in about 20% of the first 
29 frequency-transformed Galileo orbits. Galileo's or- 
bit E4 was the first orbit for which a hint of the amplitude 
modulation appeared, then in orbit G8, the signature was 
unmistakable. Other Galileo spacecraft orbits for which 
one can clearly see the modulation products are: CIO, 
El 8, and G29. Peri odograms of each of the individual 
orbits can be seen in 3raps (2001). 



1 . A strong peak near the origin, 

2. An asymmetric peak: maximum at 0.6 day -1 , center 
at 0.7 ±0.2 day -1 , 

3. An asymmetric peak: center at 1.7 ± 0.2 day -1 , 

4. A tall peak: center at 2.4 ± 0.1 day -1 , 

5. A peak: center at 3.1 ±0.1 day -1 , 

6. Harmonics of the previous three peaks, and 

7. Progressively smaller and less-defined peaks. 



2.2. Amplitude Modulation in Frequency Space 

Frequencies in frequency space can interact in numer- 
ous ways. We interpret the frequency peaks, seen in 
Fig. [H to be the result of Io's frequency of orbital rota- 
tion, Jupiter's magnetic field frequency of rotation, and 
an interaction between these two frequencies called am- 
plitude modulation (AM). The simplest case of AM is 
a sinusoid modulating the amplitude of a carrier signal, 
which is itself a sinusoid. Then the carrier signal is bro- 
ken down in frequency space into several sinusoidal os- 
cillations: x w sinfwoO + sin(woi) sin(Ot), which can 
be converted to sums of frequencies using a trigonomet- 
ric identity for sine products. The result is a signature 
in frequency space that displays a carrier frequency: cjo 
with side frequencies ("modulation products"): (ujq + £1) 
and (luq — fl). 

The process of amplitude modulation applies to Io's or- 
bital and Jupiter's rotational frequencies in the following 
way. Jupiter's rotation period is 9.8 hours corresponding 
to a frequency of 2.4 rotations per day. Io's orbital pe- 
riod (and rotation period) is 1.8 days, corresponding to 
a frequency of 0.6 rotations per day. If the dust origi- 
nates from Io, and the dust flux is modulated by Jupiter, 



2.3. Synthesizing the Frequency-transformed Dust 
Stream Data 



The main physical processes behind the frequency- 
transformed data can become clearer when one synthe- 
sizes data and compares the synthesized data to the real 
data. We have synthesized a spectrum, using typical pe- 
riods from Galileo's dust impact rate data, a Jupiter ro- 
tation and half rotation period, and an Io orbital rota- 
tional period, which show in the real data. We added 
some Gaussian-distributed noise, and after transforming 
the synthesized data into frequency space with an FFT, 
frequency peaks appear with their modulation products, 
which are in the same locations as those in the 1996-1997 
Galileo dust detector periodogram. Figures of the synthe- 
sized time series and frequency-transformed data can be 
seen in Fig. |2| 



2.4. Spacecraft Effects in Frequency Space 

Several Galileo spacecraft orbital characteristics can be 
identified in the frequency-transformed data. The first ef- 
fect is at the origin. For each orbit, the dust instrument 
receives more dust impacts while in the inner Jovian sys- 
tem than while in the outer Jovian system, which in fre- 
quency space, results in a peak at the origin. The sec- 
ond spacecraft effect is a Doppler effect between Galileo 
and Io, as the "observer" and the "source." In frequency 
space, the result is that both the Io and the Jupiter peaks 
can be smeared by Doppler shifts. The Io frequency peak 
Doppler shift is to shorter periods. This asymmetry ap- 
pears in the modulation products, as well. A tab le of the 
Doppler shift trends can be seen in praps (2001 ). 



3 



2.5. August-September 2000 Dust Storm 

In August and September 2000, both Cassini (travelling 
by at ~1 A.U. from Jupiter), and Galileo (in an orbit car- 
ried to ^250 Rj from Jupiter, where Rj=7.134 10 9 cm 
is the equatorial radius of Jupiter) detected increases 
in the rate of dust impacts, which were approximately 
100 times their nominal impact rates. In the frequency- 
transformed data from both spacecraft, Io's frequency 
signature swamped all other frequency signatures in the 
Galileo data, which was noteworthy, because both space- 
craft were located far from Jupiter, outside of Jupiter's 
magnetosphere. Galileo detected the 'dust storm' ear- 
lier during the perijove portion of its G28 orbit, during 
Days ~2 18-240. After Day 240, Galileo's impact rate 
decreased, however, Cassini observed high impact rates 
particularly on ~Day 251 and ^Day 266. (See the peri- 
odogram of the dust impact rate in Fig. |]). Cassini's dust 
detector's observational geometry is very different from 
Galileo's, however, by accident, Cassini captured the dust 
storm approximately 1 -2 weeks after Galileo detected the 
dust storm. 



2.6. Conclusions of Io's Frequency Fingerprint 

Frequency analysis via Fourier techniques of the Galileo 
dust data provides our first direct evidence of an Io dust 
source. The presence of Io's rotation frequency argues 
that Io is a localized source of charged dust particles 
because charged dust from diffuse sources would cou- 
ple to Jupiter's magnetic field and appear in frequency 
space with Jupiter's rotation frequency and its harmon- 
ics. A confirmation of Io's role as a localized charged 
dust source arises through the modulation effects. 

The Galileo dust detector periodogram data shows vari- 
abilities, orbit-to-orbit, even for some orbits which share 
similar orbital geometry. This orbit-to-orbit variability is 
a clue that, either the intervening medium or Io itself, i.e. 
its volcanoes, is a source of the variability. 



3 . JOVIAN DUST STREAM DYNAMICS 



To describe the dynamics of the dust streams, we ap- 
ply a detailed Jovian particles and fields model to sim- 
ulate a dust stream particle's trajectory as the particle 
moves from Io's orbit through Jupiter's magnetosphere 
and beyond. For the model, one needs to assume approx- 
imations for the following: 1) Jupiter's magnetic field, 
2) Jupiter's plasma, 3) dust particle density, 4) dust parti- 
cle optical property, 5) charging processes, and 6) forces. 



3.1. Jovian Dust Stream Model Details 



internal field. Addition ally, we hinge Connerne y's cur- 
rent sheet described in Connerney et al. (1981 ), to ap- 
proximate Jupiter's magnetodisk. Connerney's sheet im- 
plementation considers the magnetodisk as a perturbation 
to Jupiter's internal field. 

Jupiter's plasma is approximated using a plasma model 
which is a fit to the Vo yager 1 and 2 c old plasma mea- 
surements, described in Bagenal (1989). We assumed a 
constant mixing ratio of 50% between single ionized oxy- 
gen and sulfur ions. 

For the dust particle density, we use density values 1 .35- 
2.0 g/cm 3 . The dust particle's optical properties are man- 
ifested via Q pr , which affects the particle's dynamics 
through the radiation pressure force. For this work, a Q pr 
value for the dust particle is calculated based on the parti- 
cle's si ze, and following t he curve in Burns et al.'s classic 
paper: Burns et al. (1979] ). 



The charge of the dust particles is approximated by sum- 
ming over the currents: photoelectron emission, ion and 
electron collection, and secondary electron emission. The 
charge of the particle, which varies in time, is integrated 
simultaneously with the particle's acceleration. Here, 
we model the acceleration by considering the following 
forces: Jupiter's gravitational force, the light pressure 
force, the Lorentz fo rce, and the solar gravitational force 
(Horanyi et al. 1997). We have neglected the neutral gas 
and plasma (Coulomb) drag forces on the dust particle be- 
cause the time-scales used in these runs are short (hours 
to a day), compared to the time-scales over which those 
forces have an effect. 



More details of the Jovian dust stream model can be 
foun d in M. Horanyi's papers, such as Horanyi et al. 
(1997), and in |Grap S (2001| )1 



3.2. Dynamics Results 

A window of particle sizes exists for which dust particles 
can escape from traveling in Keplerian orbits in Jupiter's 
magnetosphere. For small dust particles, their motion is 
as plasma ions and electrons, which gyrate about Jupiter's 
magnetic field lines. For large dust particles, their mo- 
tion is governed by gravity. From Io's orbital location, 
the window of dust particle sizes for escaping particles is 
approximately radius 5 nm to 35 nm. This particle size 
is strongly dependent on the charging assumptions, espe- 
cially the secondary electron emission material assump- 
tion. In our numerical experiments thus far, we found 
that smaller dust particles can be ejected when the im- 
pacting ion or electron energy is lower than the energy of 
other impacting ions or electrons in the secondary elec- 
tron emission process. For example, a particle with sulfur 
properties can be ejected from Io's location with a min- 
imum size 4 nm, versus a particle with silicate material 
properties, which can be ejected from Io's location with 
a minimum size 6 nm. 



We approximate Jupiter's magnetic fi eld by implement- 
ing C onnerney's O4 or 06 model (Connerney 1981 



1993), which is a quadrupole expansion of the planet's 



As the dust stream particle moves through the Jovian 
magnetosphere, equilibrium potential is rarely reached. 
Therefore, as the particle moves outward, it continues 
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to collect charges, which further accelerate the particle. 
We now think that the Jovian dust stream particles move 
faster (at least 400 km/sec) t han previously assu med in 
the earlier work presented by Zook et al. (1996 ), which 
suggested dust speeds, at least 200 km/sec. 



3.3. Joint Measurements: December 2000 from 
Galileo & Cassini 



On December 30, 2000, the Cassini spacecraft closely 
flew by Jupiter, providing a simultaneous two-spacecraft 
measurement (Cassini-Galileo) of particles from one col- 
limated stream from the Jovian dust streams. Particles 
in a stream were detected by Galileo, as the spacecraft 
was orbiting inside of the Jovian magnetosphere close 
to Ganymede (8-12 Jovian radii), and then particles in 
the stream traveled to Cassini, as Cassini flew by Jupiter 
at approximately 140 Jovian radii. Figure Q shows the 
dust impact rate data for the dual dust stream measure- 
ments, the gold line denotes the Galileo rates and the 
green line indicates the Cassini rates. We assumed that 
the same dust stream at each spacecraft began where the 
black horizontal line marks the midpoint of the peak rise 
in impact rate. The travel time between the two black- 
marked peaks is approximately 7 hours. One goal of the 
dynamical modeling was to match this travel time. From 
preliminary modeling, Fig. || shows the result of one pos- 
sible trajectory of a Jovian dust streams particle released 
near Io's orbit. Here, the smallest dust particles could 
have the following range of parameters: size: 6 nanome- 
ters, density: 1.35-1.75 g/cm 3 , initial charge potential: 
1-4 V, secondary electron emission yield: 3.0, depen- 
dent on a maximum electron energy 300 eV, and a pho- 
toelectron emission yield: 0.1-1.0, which produce dust 
particle speeds: 220\450kms~ 1 (Galileo\Cassini) and 
charge potentials: 5.5\6.3 V (Galileo\Cassini). Smaller 
and larger particles than 6 nanometers result in the wrong 
direction towards Cassini, and with travel times that are 
either too fast or too slow. 
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4. SYNOPSIS 



Our work from frequency analysis shows Io as the dom- 
inant source of the Jovian dust streams. The variability 
seen in the frequency analysis shows that we might be 
able to use dust stream measurements to monitor Io's vol- 
canoes' plume activity. 

Our charging and dynamics modeling, (more details in 



jraps (2001)), shows that the dust streams' equilibrium 
potential is rarely reached in the Jovian magnetosphere, 
and that the Jovian dust stream partic les travel faster than 



found previously in Zooket al. (1996). In our preliminary 



analysis and modeling of the Galileo-Cassini dust stream 
measurements, we show that one set of conditions, which 
can match the travel times, gives a dust streams speed at 
Galileo of 220 kms -1 , with a charge potential: 5.5 V, 
and a dust streams speed at Cassini of 450 kms -1 , with 
a charge potential: 6.3 V. 
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Figure 1. A Lomb-Scargle periodogram for the first two years, 1996-1997, of Galileo dust impact rate data, with the 
specific signatures described in the text. Io 's frequency of orbital rotation, in particular, can be seen here in frequency 
space, which provides our first direct evidence that Io is a source of the Jovian dust streams. A confirmation oflo's role 
as a localized charged dust source arises through the modulation effects. The vertical solid lines mark Io 's orbital and 
Jupiter's rotational periods, and the arrows point to Jupiter's modulation products with Io straddling Jupiter's frequency. 
The first harmonic of Jupiter's rotation frequency is visible at (u)\ = 4.8 day^ 1 ) and Jupiter's modulation products with 
Io, which are straddling that first harmonic peak, can be seen, as well. The strong frequency peak near the origin at 1 
over Galileo 's orbital period is due to the Galileo spacecraft orbital geometry. 
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Synthesize^ Galileo DPS Data | plus Noise 



Noise Distribution: Goussion 
Noise Amplitude: 0.0100000 
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Figure 2. Top: A synthetic impact rate dataset, with added Gaussian noise, which includes the following periods: 48 days 
for Galileo 's orbit, 9 and 17 days for smaller data blocks, a Jupiter rotation and half rotation period, and an Io rotation 
period. Bottom: An FFT of the synthetic rate dataset. Galileo's motion (at the origin), Io's frequency of orbital rotation 
(0.6 d^ 1 ), Jupiter's full frequency (2.4 d~ x ) and half -frequency (4.8 d~ x ) of rotation can be seen, as well as the Io-Jupiter 
amplitude modulation products (sidelobes around Jupiter's frequencies). 
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Figure 3. A Lomb-Scargle periodogram of dust impact rate data from Galileo (top) and Cassini (bottom) during selected 
days in July and August 2000. The Galileo data occurred during Galileo 's G28 orbit, when the spacecraft was located 
outside of Jupiter's magnetosphere at ^250 Rjfwhere Rj=7.134 10 9 cm is the equatorial radius of Jupiter). This lo 
dust streams storm detected in interplanetary space was stronger in the earlier portion of this orbit close to perijove on 
Days 218-240 than in the later days, when Cassini detected the storm. Cassini observed high impact rates particularly 
on ^Day 251 and ~Day 266, when the spacecraft was traveling far outside of Jupiter's magnetosphere at ~7 A. U. 
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Figure 4. Dust impact rates from the Cassini-Galileo dual dust stream measurements performed on December 30, 2000. 
The marked travel time, 7 hours, is the time for a dust particle in a collimated dust stream to travel from Galileo to 
Cassini. 




Figure 5. Dust particle trajectories in x-y (top) and x-z (bottom) Jupiter-centered coordinates for different sized particles 
released from a near-Io Keplerian orbit at 6.2 Jovian radii. The colors indicate the values of the particle's charge 
potentials. The trajectory with the arrow indicates a reasonable trajectory which would bring the particle near to the 
Cassini spacecraft. 



